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We perform a precise analytic test of the instanton approximation by comparing the exact band 
spectrum of the periodic Lame potential to the tight-binding, instanton and WKB approximations. 
The instanton result gives the correct leading behavior in the semiclassical limit, while the tight- 
binding approximation does even better. WKB is off by an overall factor of -v/e/Tr. 
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Periodic quantum systems arise in many areas of physics, from crystal structures in solid state physics, to optical 
lattices in AMO physics, to solitons in polymer physics, and to the vacuum structure of QCD. An important general 
feature of these systems is the phenomenon of quantum tunneling that broadens discrete energy spectra into bands. 
Various semiclassical techniques have been developed to analyze the spectra of such systems: the tight-binding, 
WKB, and instanton approximations. In quantum field theory, the instanton approach provides important insights 
into symmetry breaking and the QCD vacuum W. In this paper we compare these approximations for an exactly 
solvable periodic system - the Lame model - with particular emphasis on the instanton approximation for the width 
of the lowest energy band. Two well-studied quantum mechanical instanton models |0-Hl are the double-well, V{(1>) 



i0^(l — y/g(t))'^, and the 'Sine-Gordon', V{(f>) = -{1 — cos(y^0)). The Lame model, V{(f)) = ^sn^(^0|i/), provides 



a new example, and has the advantage that it can be solved exactly without the instanton approximation, while the 
instanton calculation can also be done analytically. We find that there is a difference between the semiclassical and 
tight-binding limits, limits that are often regarded as synonymous. We also highlight a new connection between the 
non-perturbative instanton approach and the algebraic approach to spectra [^,|6| . 
Consider the Schrodinger-like Lame equation 






l)i/sn2(0|i/)*(0) =£*((/.) 



(1) 



where j = 1, 2, 3, . . . is an integer, sn(0|i/) is one of the Jacobi elliptic functions, and < i^ < 1 is the elliptic parameter. 
This equation is known to be exactly solvable in terms of theta functions, and is known to have a spectrum of j bound 
bands 0. The function sn^{(j)\v) is periodic, with period 2K{h'), where K{i/) = J^ d9/\/l — i/sin^ 9 is the elliptic 
quarter period. For small i/ (say v < 0.2), sn^((/)|;/) looks to the eye like sin^(0). For large f (say v > 0.9), sn^((/)|z^) 
looks like a sequence of periodically displaced Poschl- Teller — sech^(0) potentials (see Fig. 1), displaced with a period 
that diverges logarithmically, 2if (i/) ~ log( j^), as z^ ^ 1. 




FIG. 1. The function sn^ 



with u = 0.999, as a function of (f). The period is 2ii:(0.999) = 9.68. 



While the Lame equation is exactly solvable for any j and ly, explicit expressions for the wavefunctions and energies 
are cumbersome for j > 2 . A dramatic simplification is that the band-edge energies of (|l|) are given by the 2j -I- 1 
eigenvalues of the finite dimensional matrix Jj -I- vJy, where Jx and Jy are su{2) generators in a spin j representation 
|g|,g| . Thus we can evaluate the width AS of the lowest band as the difference between the two smallest eigenvalues 
of the matrix J^ -|- vj'^. For example, for j = 1, Af <=''^'=* = 1 - J^; for j = 2, Af '=''^'=* = _l _ jy -^ 2^^2 
j = 3, AS''''^''^ = 3(1 -1^) + 2Vl - i^ + 4i/2 - 2\/4~^^^ 



u; for j = 2, Af '=''^'=* = _1 _ jy -^ 2Vv^ - i/ + 1; and for 
2. We will be interested later in large values of j, in which 



case such explicit expressions, as functions of v, become more difficult to derive. Instead, we have shown algebraically 
that for any j, the exact energy splitting of the lowest band, as i^ — > 1, is 
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8jr(j + l/2) 

4-'"V^r(j) 



{1-i.y {i + l-l{i-^) + 



(2) 



As 1/ ^ 1 the period becomes infinite, which suppresses tunnehng, and we therefore expect this to be relevant for 
the tight-binding and semiclassical approximations. Note that the band- width (0) vanishes as z^ ^ 1, and vanishes 
more rapidly for larger values of j. We now compare this exact result (||) with the tight-binding, instanton and WKB 
approximations. 

In the tight-binding approximation of solid state physics [BJ the separation L between neighboring wells becomes 
large and we treat the potential as a sum of periodically displaced 'atomic' potential wells: F(0) = ^„ [/((/) — nL). 
In the Lame case, this physical approximation is explicitly realized by the remarkable mathematical identity pi| : 

As expected, each 'atomic well' has the form of a Poschl- Teller well (but the rescaling factor i^ is non-obvious). 
Including the jij + 1) factor from (||), each atomic well has j discrete bound states, and the effect of the periodic sum 
is to broaden these states into the j bound bands of the Lame potential. For the lowest band we use the ground state 

^o{<j)) = V 2K''r(/) ^sch-' (^^(/)) of the atomic well. The width of this band can be calculated using standard solid 
state techniques 

/C30 
dct> J2 sech^ {^i^ - 2"^)) *o(0)*o(0 - 2K) 

where in the second line we have kept dominant terms as ;/ — > 1, and used the fact that exTp[—TrK{i')/K'{iy)] ^ 
"1^(1 + 5(1 — J^) + • . .). This result (Q) agrees precisely with the exact result (||) to this order in 1 — i^. Therefore, for 
any j, the tight-binding approximation is good as i^ ^ 1; i.e. as the separation between atomic wells becomes large. 
We now turn to an instanton evaluation of the width of the lowest band, expected to be good in the semiclassical 
limit in which tunneling effects are small. Naively, one might expect that this is also just the v ^ I limit in which the 
wells become infinitely separated, but it is actually more interesting than this. To make contact with the standard 
instanton approach [^pQ we define 

V{<f>)^^^sn'{^cj^W)^^^'^g^^^cP^ + ... (5) 

where g is some coupling (which we will relate to j and v below) , and we have chosen units so that the perturbative 
mass in a given well is 1. The width of the lowest band can be found by considering the Euclidean path integral 
connecting two neighboring classical minima of the potential; here </) = and (jj = 2K/y/g for example. Rescaling the 
field variable to x = V5'/', 



exp^--5[0]j =exp^--y dt [-i^Y + -sn^ M.)] ) (6) 

Thus the semiclassical limit is hg -^ 1, or -^h -^ l^ioak '■ the ground state in each well is far below the barrier height. 
There is a standard technique for computing the instanton approximation for the width of the lowest band 1^-^] . In 
the semiclassical limit, the Euclidean path integral is dominated by 'instanton' solutions to the Euclidean equations 
of motion, satisfying the first-order equation Xinst = ^2V^(xinst)- Using the rescaled potential ^(x) — ^^.u^ixW), we 
find the Lame instanton: 

Xi„st(i) = K{iy) + sn-i(tanh(i)|i.) (7) 

Here the integration constant has been chosen so that the instanton is centered at i = 0, and sn^^ means the inverse 
function (a standard function in Mathematica). This instanton interpolates between x = at i = — cxd, and x = '2K{i') 
at i = -|-oo, as shown in Fig. 2. The corresponding Euclidean action is 



So^J dxv/27M=-=logf- 



(8) 




FIG. 2. The Lame instanton Xi.ist(t) in (0), plotted for v = 0.9. Note that 2K{0.9) = 5.15618. 

The leading exponential factor in the instanton expression for the band width is exp[— 5'o/(?i5)]. But there is also a 
prefactor that is related to the determinant of the fluctuation operator — -^ + V" {xinst{t)) ■ This prefactor is physically 
significant, as it encapsulates the collective coordinate effects of fluctuations about the instantons. Here, 

Vtanh*(t) + 2{l-iy) tanh^(i) - l" 



V'Uinstit)) ^(l-i^) 



(l-i/tanh^(i))2 



(9) 



which is plotted in Fig. 3. We see that V" {xinst{t)) is highly localized in the vicinity of the instanton, and tends to 1 
(the square of the mass) as t —^ ±oo, as expected. 




FIG. 3. The fluctuation potential V"(xinst(i)) in ^, plotted for u = 0.9. 

Coleman has given a very simple method y] for evaluating the fluctuation determinant (with zero mode removed) , 
resulting in the formula (for a periodic potential) 



y'Trhg rig 



(10) 



where the constant c is simply determined by the asymptotic behavior of the zero mode: Xinst(0 '^^'ce l*las|i|— s-cx). 

,,^ e~l*L Thus we obtain the Lame instanton result 



Here, the zero mode is Xinst(0 = sech(t)/y 1 — z/tanh [t) 
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(11) 



To relate this to the eigenvalues £ of the Lame equation dl^) , we compare (nf) with the Schrodinger equation, 
^sn2(^0)* = E'^, for the potential V{<j)) in (||). We therefore identify 



2 



^j(i + i) 



1 



f = 



2E 



Thus the instanton approximation for the band-width of the eigenvalue £ in (fil) is 



r* + 



(12) 



As a function of j this differs from the exact and tight-binding expressions (0) and (H) , even in the large period limit 
V ^ \. However, from (ph and (O), the semiclassical limit hg <^ 1 means iyj{j -f 1) ^ 1. So, to compare the instanton 



formula ( p^ ) with the results (g) and (Q) it is not enough to take z/ ^ 1; we also need to take j to be large. Physically, 
it is not enough to take far-separated wells; they must also be deep wells. For large j and i/ near 1, 

Ar-* ~ -^ (1 - i^y (^1 + 1^(1 - ^) + . . .) (14) 

which agrees perfectly with the large j limit (using Stirling's formula) of the exact result (^. 

As another test of our instanton formula ( |l3|) we can fix i' to any value (not necessarily near 1) and take j large 
(in order to be in the semiclassical regime), and compare to the band- width obtained numerically from the two 
lowest eigenvalues of the matrix J^ + i^J^. The results are shown in Fig. 4, showing 10% agreement for j > 10. In 
terms of the algebraic spectral program [plpj , we find that non-perturbative instantons play an interesting role in the 
semiclassical limit via eigenvalue differences for finite dimensional, but very large matrices. As a final calculation, 
we have computed the WKB approximation for the energy splitting, and we find that in the semiclassical limit 
^^WKB _ 2^ gxp[_ J_ J_^p dx\/sn^{x) ~ %] = x/^^E™^*-, confirming that WKB gets the correct leading exponential 
but the WKB prefactor normalization is incorrect |12|. 



FIG. 4. The ratio of the exact numerical band-width to the instanton expression (|l3|), for various values of j, with u = 0.5. 

To conclude, the instanton approximation gives the correct leading semiclassical result for the lowest band-width. 
But the Lame model (Q) has two independent parameters, i' and j, that allow us to probe separately the 'far-separated- 
weir and 'deep-well' limits respectively. The tight-binding approximation (something of a misnomer) is good for large 
period, for any well-depth; while the instanton approximation requires deep wells for a given period so that the 
combination h'j{j + 1) is large. It would be interesting to study instanton - anti-instanton interactions and correlation 
functions in the Lame system, and to compare to lattice simulations in the spirit of |T3]. Also, Lame solitons with 
spatial profile of the form (uh will arise in the corresponding 1 + 1 dimensional model. 
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